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ABSTRACT 



The spatial distribution of 3.3 |J.m PAH and associated emission in the 3.3" x 6.0" 
inner region of the Red Rectangle nebula has been determined using the UIST imager- 
spectromcter at the United Kingdom Infrared Telescope (UKIRT). Interpretation of 
the 3.3 |J.m feature as comprising two spectroscopic components centred at 3.30 |i.m 
and 3.28 |xm, as put forward by Song et al. (2003, 2007), is supported by these data 
which reveal that they have different spatial distributions. It is deduced that there are 
two classes of 3.3 |a.m band carrier with a peak wavelength separation of ~0.02 |j.m. 
From comparison of the 3.3 |J.m observations with laboratory and theoretical spectra 
for a range of PAH molecules it is proposed that the 3.28 |J.m and 3.30 |J.m components 
arise from 'bay' and 'non-bay' hydrogen sites, respectively, on the periphery of small 
neutral PAHs. Observational data are also obtained for i^band continuum emission 
and for the Pfund e hydrogen recombination line. 

Key words: stars: individual (Red Rectangle) - techniques: spectroscopic- ISM: 
molecules - ISM: lines and bands - ISM: abundances 



1 INTRODUCTION 

The 3.3 [im emission band is a well-known member of the 
family of unidentified infrared 'UIR' or 'Aromatic IR' bands 
(AIBs) that are attributed to polycyclic aromatic hydro- 
carbon (PAH) molecules. It is the shortest wavelength UIR 
band and is assigned to the infrared active C-H stretch- 
ing modes. This feature, in common with other major UIR 
bands at 6.2, 7.7, 8.6 and 11.2 [i-m, is ascribed to PAHs 
at a generic rather tha n spec ifi c molecu l ar lev el - for dis- 
cussio n see: iTokunagal (| 19971 ): iTielend (|2008D : IWalmsIevI 
ters ( 20111). Only in the ca ses of benzene (CeHe) 



(|2009h : |Peete 



in CRL 618 JCernicharo et all [iooH ) and SMP LMC 11 
iBernard-Salas et all 120061') . Ceo and C70 in the planetary 
nebula Tc 1 (|Cami et al.ll2010l ') a nd Ceo in the reflect ion neb- 
ulae NGC 7023 and NGC 2023 (|Sellgren et al.ll2010l ) and in 
other objects, have specific Ce - ring molecules been identi- 
fied. 

The biconical Red Rectangle nebula, centred on the bi- 
nary star HD 44179, is a remarkable object with a very in- 
tense display of UIR bands. Long-slit studies of the 3.3 \im 
emission feature of the Red Rectangle and other objects 



(|Song et al.ll2003l . [2007l '). have suggested that the band pro- 
file can be interpreted in terms of two emission compo- 
nents centred at 3.28 (im and 3.30 (i-rcQ. Exposures with a slit 
aligned along the NW whisker of the nebula revealed that 
3.30 [im emission appears on-star and that the 3.3 \un fea- 
ture changes shape with offset due to addition of a 3.28 \im 
component. Although this two-component description works 
well as a basis for interpreting the spectra, it does not pre- 
clude weak emission at slightly different wavelengths and 
there are indi cations of some w avelength variation as a func- 
tion of offset (|Song et al.ll2007l ). 

In this Integral Field Unit (IFU)-based study the spa- 
tial distribution of the 3.3 \im band is investigated over the 
inner part of the nebula (3.3" x 6.0"), centred on the star 
HD 44179. The aim is to improve knowledge of the distri- 
bution of PAHs in the Red Rectangle and of the size, shape 
and charge state of 3.3 \im band carriers both in this object 
and more generally in Galactic and extragalactic sources. 
The paper is arranged first with a review of previous ob- 
servational, experimental, theoretical and modelling studies 
of the 3.3 [xm band (Section 2), followed by observational 
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^ As in I Song et al. I 1I2OO7I ') we use the nominal description '3.3' 
to refer to the whole 3.3 \im profile, with the wavelengths of the 
two components derived from profile fitting written with a second 
decimal place viz 3.28 and 3.30 |J.m. 
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details (Section 3) and the derived IFU images and their 
spectral and spatial decomposition (Section 4). In Section 
5 we consider the origin of the 3.3 ^.m band with particular 
reference to experimental gas-phase spectroscopic data for 
small PAHs and their bay/non-bay character, followed by 
DFT calculations of vibrational frequencies and transition 
strengths and the incorporation of this experimental and 
theoretical information into a UV-pumped emission model 
(Section 6). Comparison is made between the 3.3 [im astro- 
nomical images and those for 'blue luminescence' (Section 
7), and a summary of results is given in Section 8. 



2 OVERVIEW 

In order to provide a backdrop for the following Sections 
the current state of knowledge of the 3.3 (xm PAH feature is 
summarized here from observational, laboratory, theoretical 
and modelling perspectives. 



2.1 Astronomical studies of the 3.3 vim band 

Using ground-based studies iTokunaga et al.l (|l99ll ) made 
the first classification of the 3.3 ^m band and identified a 
common Type 1 and a rare Type 2. Following ISO obser- 
vations, an (A, B) classification scheme based on the band 
shape and central wa velength was put forward for a large 
number of objects by Ivan Diedenhoven et al. I (|2004 ). with 
designations A3. 3, BI3.3 and B23,3. The A3. 3 set is com- 
mon, has an approximately symmetrical profile with peak 
position of - 3.290 ^m and a FWHM of 0.040 ^m; it is 
similar to Type 1. In their two-component interpretation, 
ISong et all (|2007l ) showed that this corresponds to approxi- 
mately equal contributions of 3.28 [im and 3.30 [im 'carriers', 
which results in a peak wavelength of ~ 3.29 \im. 

A smaller group of objects was classed as BI3.3 
with maximum intensity at slightly longer wavelength ('^ 
3.293 ^im) and with a somewhat narrower width than for 
A3. 3. In the Song et al. (2007) interpretation this corre- 
sponds to a smaller relative contribution from the 3.28 \im 
component. The class B23.3 has only one member and is 
found for ISO spectra of the Red Rectangle. However, these 
data comprise a 14" x 20" exposure of both star and nebula 
and the profile in fact varies within the object (|Song et al.l 
l2003f ): the peak wavelength is yet longer at ~ 3.297 [im 
(FWHM ~ 0.037 \im) and lies fairly close to a 'p ure' 3.30 [im 



compo nent as found for on-star PAH emissi on bvlSong et al 



2003 ) . It is similar to the rarer Type 2 of iTokunaga et al 



19911 ). The overall 3.3 vim profile in the Red Rectangle de- 
velops from Type 2 (~ B23.3) on-star towards Type 1 (or ap- 
proximately) BI3.3 at ~ 7.5" ofi'set. Expressed in terms of 
spectral components, for the Red Rectangle the Lorentzian- 
shaped 3.30 vim feature dominates on-star and the 3.28 vim 
feature grows (relatively) with ofi'set, reaching about 50% of 
the 3.30 vim intensity at ~ 7.5". 

The general sequence A3. 3 — BI3.3 — B23.3 can be 
summarized semi-quantitatively in terms of a 3.28/3.30 in- 
tensity ratio as ^ 1, ~ 0.5 and ^0.2, respectively, although 
a continuous range of ratios is very likely to exist among dif- 
ferent astronomical sources, and even within a single object 
as occurs in the Red Rectangle. 



2.2 Laboratory studies of the 3.3 vim band 

Numerous experiments have been conducted in absorption 
on individual PAH molecules trapped i n cold inert gas matri- 
ces, e. g. neutral PAHs w ith 2-4 rings pludgins fc SandfordI 



Il998al ) and > 5 rings jHudgins fc SandfordI Il998bh . Al- 



though somewhat removed from gas-phase PAHs under in- 
terstellar conditions because of the extremely low temper- 
ature and matrix perturbations, these measurements al- 
low easy comparison of infrared spectra of a wide range 
of both neutral and ionised PAHs. The most extensive ex- 
perimental information for the 3.3 vim band of PAHs is 
iven in the NASA Ame s Infrared Spectroscopic Database 
Bauschlicher et al.llioiol '). 

Experiments have been conducted on the 3.3 vim 
vib rational band of gas-phase PAHs i n both absorption 
e.g.lSemmler. Yang fc Crawford (liggil): IJoblin et all (Il994l . 
■ N 19961. 11997(1 



1995h: ICaiie et al 



Pirali et al.1 
20 111) 



and 



( 20091 ') 

in emission 



\): iHunevcutt et al.1 (l2004h: 



NIST MS Data Center / S E SteinI 



Schlemmer et al.l 



Williams fc Leond |[l995l ) : ICook et al.l (Il996l. 



Wagner. Kim fc Savkallvl I2OO0I I: iKinTfc Savkallvl (|200a ). 
In a close simulation of interstellar conditions, UV laser 
excitation of cold gas-p hase naphthalene in a supersonic jet 
(|Schlemmer et alii 19941 ') yielded a single-peaked IR emission 
band centred at 3.30 vim (see their Figure 1 A); notably 
it is significantly red-shifted with respect to the gas-phase 
absorption spectrum at ^500 K which p eaks at 3.26 vim 
INIST MS Data Center / S. E. SteinI [20l3) and a jet-cold 
gas-phase spectrum which sh ows two reso l ved pe aks at 3.26 
and 3.25 vim. As shown by Ijoblin et all (|l995l ). it is not 
expected that gas-phase C-H absorption wavelengths at 
low temperatures will match astronomical observations as 
there is a marked dependence of the peak wavelength on 
temperature. Given that the astronomical infrared emission 
is necessarily from significantly heated PAH molecules 
following UV/visible photoexcitation, in this work we have 
chosen to compare astronomical spectra with those of 
heated PAHs, where available, and this is discussed further 
in Section 5. 

Experimental work on condensed-phase carbon mate- 
rial in both absorption and emission reveals a C-H stretch- 
ing feature which matches quite wel l the astrophysical data 
JSakata et al.lll9 90': 'Goto ct al."l997'; 'Scott. Dule v fc Jahanil 
ll997l : lHu fcplll cv 2008; Dulcy & Wilhams 2011). However, 
the physical state of the material is significantly different 
from the gas-phase PAHs which are more commonly invoked 
as UIR band carriers. 



2.3 Theoretical studies of the 3.3 vim band 

Although high-level ab-initio calculations of C-H vi- 
bratio nal frequencies and intensities have been under- 
taken ("de Frees et al.|[]~993l : iPauzat. Talbi fc Emngeijll997l : 



IPauzat 2011 ). these are necessarily limited to small sys- 
tems. The overwhelming majority of computations have 
been made using de nsity functional theory (DFT) mostly 
with B3LYP/4-31G (lLanghofj|l996l : iBauschlicher fc Bakel 
I2OOOI : iBauschlicheil I2OO2I ). The most extensive theoreti- 
cal information for vibrational transitions for PAHs is 
given in the NASA Ames Infrared Spectroscopic Database 
jBauschlicher et al.1 [ioiol ). A complementary database by 



© 2012 RAS, MNRAS OOO.HlfTOl 



iMalloci. Joblin fc Mul3 (|2007l ') includes energetic, rota- 
tional, vibrational and electronic properties of PAHs in a 
wide range of charge states. 

In a recent study, vibrational frequencies for large 'reg- 
ular' and 'irregular' PAHs have been calculated for PAHs 
with up to 130 carbon atoms includ i ng th e 3.3 [i m band 
iBauschlicher. Peeters fc Allamand"oial |2008| . |2009| ). These 
authors conclude that cationic, anionic and neutral PAHs 
all contribute to the 3.3 [im feature, the anion to neutral ra- 
tio being estimated to be ~ 0.7. It was suggested that the 
relative contributions of anionic and neutral forms for com- 
pact PAHs may be responsible f o r the Type 1 and Type 2 
band profiles of iTokunaga et abl (|l99lh . where the shorter 
wavelength contribution would be due to neutrals. It was 
also concluded that 'the vast majority of astronomical PAHs 
responsible for the emission in the 3 [im region are compact, 
not irregular, elongated, or bent'. This deduction was based 
on the relative flux of the continuum-subtracted 3.3 [un 
feature to the residual above the continuum at 3.23 and 
3.24 \im, respectively, for neutrals and anions, and where a 
15 cm~^ red shift was invoked. 
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Table 1. Wavelength ranges selected for the spectral images. 



Feature 


Origin 


Range (|xm) 


L-band continuum 


Dust 


z.5oU-3.z4U; 
3.540-3.610 


3.3 [xm (total) 


C-H stretch 


3.260-3.350 


3.28 & 3.30 iJ-m 


C-H stretch 


Deconvolution 
(see text) 


3.4 |xm 


C-H2,3 stretch 


3.375-3.427 


Pfund e 


H recombination 


2.980-3.100 



2.4 Modelling studies of 3.3 (xm emission 

A number of studies have been directed towards under- 
standing astronomical 3.3 \im emission, u sually as one of a 
wider r a nge of UIR bands - se c for example Cook fc Savkallv 



Pech. Joblin fc Boisscl (,20o3 ): IPathak fc Rastogi 



iBauschlicher et al.l ( 2O10l ). These studies show how 



the 3.3 [im band is dominated by small PAHs, which can 
reach peak temperatures high enough to allow emission 
in hi gh frequency modes (|Schutte. Tielens fc AUamandolal 
19931). The study of most relevance to this work is that of 



Mulas et al.l (|2006l ) in which infrared emission fluxes for an- 
thracene, phenanthrene and pyrene were computed in de- 
tail with reference to the radiation field in the Red Rect- 
angle. Stimulated in significant part by the discovery of 
'Blue Luminescence' (BL) and the suggestion that this is 
due t o emission from molecules with th r ee or four aromatic 
units (|Viih. Witt fc Gordonll2004l . l2005h . [Mulas et all (|2006l ) 
showed that the computed 3.3 vim fluxes for anthracene, 
phenanthrene and pyrene were reasonably consistent with 
ISO data for 3.3 yun emission of the Red Rectangle. The 
possible link between BL and the 3.3 [im PAH band is dis- 
cussed further in Section 7. 



3 IFU OBSERVATIONAL DETAILS 

HD 44179 was observed on three separate occasions, 2005 
February 19, February 22 and March 5, using the UIST 
imager-spectrometer at the 3.8 m United Kingdom Infrared 
Telescope (UKIRT). The observations were carried out us- 
ing the integral field unit (IFU) in combination with the 
short-i grism, giving a central wavelength of 3.27 \im and 
a wavelength coverage of approximately 0.7 [im. The UIST 
IFU covers a rotatable 3.3" x 6.0" region of sky with a plate 
scale of 0.12" X 0.24". Light incident on the IFU is then fed 
to the 14 useable 0.24" x 6.0" slitlets, then onto the dis- 
persing grism (R ~ 1400) and finally onto a 1024 x 1024 
pixel InSb array. A spectral image, individual spectra and a 
"white-light" image can thus be derived from the data. 



HR 2491 (AlV, K = -1.39) and HR 1713 (B8 lab, 
K — 0.213) were observed at matching airmass both before 
and after the observations of HD 44179 in order to remove 
telluric features and for flux calibration. The image rotator 
was used to align the long-axis of the IFU with the outflow 
from HD 44179 (position angle 13'' E of N) and the tele- 
scope was nodded approximately 7' NW for sky frames due 
to HD 44179's extended nature. An ambiguity was noted 
in the recorded coordinate system which affected the N- 
S orientation. Fortunately this was completely resolved as 
the emission flux in the Southernmost lobe in two earlier 
IR imaging studies has been shown to be about twice as 
high as for the Northern equivalent ijMekarnia et al ] |l998l : 
iMen'shchikov et al.|[l99i ). On the first two nights 10 x 2-sec 
exposures were taken at each nod position which was later 
optimised to 2 x 8-secs for the final night giving a total on- 
source exposure time of 4544 s. An internal halogen black- 
body source was used to fiatfield the data and wavelength 
calibration was achieved via identification of sky lines. Data 
were reduced initially using UKIRT's ORAC-DR data re- 
duction pipeline and then further reduced and analysed us- 
ing FIGARO, KAPPA and DATACUBE routines in the STAR- 
LINK software distribution. The i-band seeing for 19 and 22 
February was approximately 0.7" and sub-0.5" for 5 March 
2010. For ease of discussion, in the following text the de- 
scription 'N' and 'S' is used even though the orientation of 
the biconical axis is tilted 13° with respect to the actual N-S 
axis. 

Table [l] gives the wavelength ranges used in produc- 
ing spectral maps for the continuum and the PAH bands 
at 3.3 vim. The spectral maps for the '3.28' and '3.30' [im 
components were produced by performing a robust double- 
Lorentzian fi t to the full 3.3 yi m band using the IDL routine 
MPFITEXPR ()Markwardtll2009l ). The 3.4 jj,m band which ap- 
pears in spectra offset from the central star was also consid- 
ered. The fiux-calibrated spectrum produced by the pipeline 
showed the presence of the Pfund e hydrogen emission line 
at 3.04 /im, for which a continuum-subtracted image was 
also produced. 
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1.5 3.0 4.5 6.0 0.6 1.2 1.8 2.4 




Continuum (L-band) 3.3 vm band 



Figure 1. Left panel. Spectral map (3.3" X 6.0") of the contin- 
uum emission, integrated in the L-band. The colour scale is in 
units of 10-13 Wm-2 

lJ.m ^. The contours represent levels from 
5 to 95% of the maximum flux, with steps of 10%. The intensity 
maximum is centred on HD 44179 (blue cross) and extends grad- 
ually in the outer region as an ellipse with major axis 1.35" and 
minor axis 1.1". Right panel. Spectral map of the 3.3 [xm band, 
with colour scale in units of 10—^^ 'V^ra~^\i.m~^ . The contours arc 
as in the left panel. The contours form a bipolar shape, with the 
Southern lobe somewhat brighter. Colour versions of all figures 
are available in the on-line paper. 



0.35 0.70 1.05 1.40 0.55 0.70 1.05 1.40 




3.28 Mm component 3.30 |jm component 



Figure 2. Left panel. Spectral map (3.3" X 6.0") of the 3.28 urn 
component of the 3.3 \aa emission of the Red Rectangle. The con- 
tours represent levels from 15% to 95% with steps of 10%. Two 
bright peaks (separated by 1.3") appear north and south of the 
central star, but there is no emission from the centre. Right panel. 
Spectral map of the 3.30 |xm component; the contours correspond 
to levels from 5% to 95% in steps of 10%. The image is simi- 
lar to that of the 3.3 \ua band. The colour scale is in units of 
10^11 Wm^^ixm^i for both maps. 



4 IFU IMAGING RESULTS 

4.1 Continuum, 3.3 (xm band and the Pfunde line 

Figure 1 shows the spatial distribution of i-band contin- 
uum emission and the continuum-subtracted 3.3 (xm (to- 
tal) feature. The continuum emission is centrally peaked on 
HD 44179 and slightly extended along the bicone ('N - S') 
axis. In contrast the overall 3.3 \im band displays a pro- 
nounced bipolar struc ture. It may be compared with that of 
iMekarnia et al.1 (| 19981 ) for the inner c. 2"x 2", noting that 
in their study the data were obtained with a filter centered 
at 3.28 n-m with a spe ctral width of 0.16 \im. The image of 
IMekarnia et al.l (| 19981 ) diflers significantly from the 3.3 ^.m 
image in Figure 1 because it also contains continuum signal 
that was not subtracted and which was estimated to con- 
tribute c. 60% of the PAH emission signal. Figure 1 shows 
that there is some 3.3 \an PAH emission on-star (within the 
seeing constraint), with extended emission peaking c. 0.7" 
along both directions of the bipolar axis. No tably the South- 
ern m axim um is stronger, as in the ima ges of lMekarnia et al] 
1 19981 ) and lMen'shchikov et all (| 19981 ). This has been inter- 
preted as being due to inclination of the entire object which 
causes the Northern lobe emissio n to be partially obsc ured 
by dust in the circumbinary disc (|Mekarnia et al.ll 19981 ). 

The results obtained from the double Lorentzian de- 
composition of the 3.3 \ira band are shown in Figure 2. These 
spatial distributions differ in the inner region with negligible 
3.28 \im - type emission directly towards HD 44179. Inten- 



X 3.3|jm 
X 3.30 comp. 
X 3.28 comp. 



S«5 



-1.0 0.0 

Offset (arcsec) 



Figure 3. Normalised intensity profiles derived from the spec- 
tral maps, cut along the major axis at P. A. = 13°. Comparison 
between the N-S profiles for the 3.3 |J.m band, and the 3.30 and 
3.28 |J.m components is shown. 



sity profiles for the 3.3 \xm feature, and the '3.28' \xm and 
'3.30' \m\ features are presented in Figure 3. The 3.28 \xm - 
type emission grows rapidly with increasing offset from the 
central star and peaks at the same positions (N and S) as 
for the 3.30 ^m emission, so apart from the innermost re- 
gion the two images are quite similar. It appears that a 
new form of PAH develops with increasing distance from 
the star, resulting in a set of PAHs that have C-H emission 
at a shorter wavelength. The difference in the distributions 
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0.4 0.8 1.2 1.6 0.95 1.90 2.85 3.80 




3.40 Mm band pfund E band 



Figure 4. Left panel. Spectral map of the 3.4 \ira emission band. 
Colour scale is in units of 10~^^ Wm~^|^m~^. The contours corre- 
spond to 3, 6, 9, 12, 22, 32, 42 and 52% of the maximum flux value. 
The 3.4 |j.m emission is not present in the vicinity of the central 
star and has a bipolar distribution. The Southern lobe appears 
stronger than the Northern equivalent. Right panel. Spectral map 
(3.3" X 6.0") of the Pfund e hydrogen emission line. The colour 
scale is in units of 10-^*^ Wm-2 |j.m . The contours correspond 
to levels 5% to 75%, in steps of 10%. 




non-bay 



Figure 5. Illustration of bay (four) and non-bay (eight) hydrogen 
sites in the molecule perylene, C20H12. 

h as four bay and eight non-bay hydrogens. The discussion 
of Ivan Diedenhoven et all (|2004l ) was concerned principally 
with the influence of bay/non-bay hydrogens on overall line 
width rather than band frequency and with specific ref- 
erence to low-temperature matrix data which have multi- 
ple components. Based on the classifications of the 3.3 vim 
and 11.2 (xm bands among different objects, it was con- 
cluded that molecular edge structure 'is likely not the main 
cause of the variation in the F WHM of the 3.3 \im feature' 
(Ivan Diedenhoven et al.|[20o3 ). We now examine 'bay' and 
'non-bay' sites as a basis for explaining the 3.28 and 3.30 \im 
spectral components with reference to laboratory and theo- 
retical data. 



seen in these IFU images provides evidence for two types 
of 3.3 (i.m carrier, and supports the earlier conclusion from 
long-slit spectrosco pic studies of the 3.3 \im band profile in 
the Red Rectangle (|Song et al.ll2003l . |2007| ). 

The image (Figure 4) for the continuum-subtracted fea- 
ture near 3.4 [im (see Table 1) is similar to that of the 
3.28 (xm compon ent and compleme nts earlier spectroscopic 
observations by ISong et al.l (|2003l ). neither feature being 
present near to the star. The 3.4 fi.ni emission is assigned to 
the C-H stretch arising through addition of a methyl group 
(CH3) to the PAH edge or by hydrogenation as to make an 
sp^ - bonded c arbon. Althou gh not completely understood 
as discussed bv iTielen^ ()2008l ). there is little doubt that the 
3.4 fim feature arises from a chemical addition to a standard 
aromatic PAH structure. Figure 4 also shows the distribu- 
tion of Pfund e hydrogen emission which traces the ionised 
central region of the nebula. 



5 BAY AND NON-BAY CHARACTER OF 
SMALL PAHS 

Pr ofile variation in the 3. 3 ytm f eature has been discussed 
by Ivan Diedenhoven et al.l (|2004l ) with reference to charge 
state, size, edge structure and heteroatom substitution. In 
particular these authors commented on differences that 
could arise according to whether the relevant hydrogen 
atoms occupied 'non-bay' sites, as in very symmetrical com- 
pact PAHs, or 'bay' sites where nearby C-H stretching mo- 
tions could interfere; an example illustrating bay and non- 
bay sites is given in Figure 5 for the molecule perylene which 



5.1 Experimental data 

5.1.1 3.3 nm PAH emission following UV-pumping 

Infrared emission from PAHs in astrophysical environments 
is from vibrationally hot molecules at an effective tempera- 
ture which depends mostly on the energy of the absorbed 
photon and molecular size. Laboratory infrared emission 
data obtained from from UV/vis pumped PAHs are optimal 
for comparison with astronomical spectra but the number 
of PAHs studied in this way is, unfortunately, small and the 
data cannot readily be used in evaluating the possible role 
of bay/non-bay character in relation to the 3.3 \im band. 



5.1.2 3.3 PAH absorption at low temperature 

Considering first experiments on PAHs in low-temperature 
inert gas matrices, the C-H stretching spectral region is 
known to be sensitive to the degree of b ay /non-b a y char - 
acter, as re marked on by Ivan D icdcnhp ven et al.l (|2004h . 
Figure 1 of iHudgins fc SandfordI (|l998ai ) and Figure 2 of 
iHudgins fc SandfordI (|l998bl ). show that molecules with non- 
bay (nb) hydrogens only, such as pyrene with 10 non-bay 
hydrogens (CieHio; lOnb), have a single narrow dominant 
absorption feature with peak maximum at ~ 3050 cm^^ 
in solid argon. For a less symmetrical molecule with some 
bay (b) hydrogens such as chrysene (C18H12; 8nb -|- 4b), 
the spectral feature is broader with the C-H feature centred 
~20 cm~^ to higher frequency. The triphenylene (6nb -I- 6b) 
and hexabenzocoronene A (C42Hig; 6nb -f 12b) molecules 
have many more bay hydrogens and their strongest absorp- 
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Table 2. Gas-phase absorption wavenumber s for PAHs from 
Table 2. of lSemmler. Yang fc Crawford! digQll ) in which relative 
band strengths are given. The (nb, b) labelling and presentation 
of the wavenumbers in two columns is our work. The molecules 



Wovelength (/um) 
3.21 3.25 3.29 



3.33 



in rows 3-6 are the same as those in Figure 6. 


Molecule 


Formula 


nb 


b 






Anthracene 


C14H10 


10 





3063 


- 


Phenanthrene 


C14H10 


8 


2 


3065 




Pyrene 


CieHio 


10 





3051'' 




Benz [a] anthracene 


C18H12 


10 


2 


3062 




Chrysene 


C18H12 


8 


4 


3066 


3089 


Triphenylene'^ 


C18H12 


6 


6 




3084 


Bcnzo [c] phenanthrene*^ 


C18H12 


10 


2 


3058 




Bcnzo [a] pyrene 


C20H12 


10 


2 


3056 




Benzo [e] pyrene 


C20H12 


8 


4 


3056 


3086 


Perylene 


C20H12 


8 


4 


3059 


3087 


Dibenz[a,c] anthracene 


C22H14 


10 


4 


3066 


3078 


Dibenz[a,h] anthracene 


C22H14 


10 


4 


3064 


3077 



Notes, a. Wave number units, b. 3052 cm at 570 K 
jjoblin et al.lll994l'l . c. A weaker band at 3032 cm ^ is reported. 
d. The interacting Hs have a different position (fjord) from the 
other cases in the Table; a weak band at 3013 cm^^ is reported. 

tion is at around 3100 cm^^ in the matrix environment. 
These spectral data indicate that the presence of 'bay' sites 
introduces higher frequency absorption. There also exist 
data for isolated jet-cold neutral gas-phase PAHs which, 
for the five PAHs in common, generally reflects the appear- 
ance of the matrix da ta with frequency agre ement usually 
to within a few cm~^ (|Hunevcutt et aLll2004h . 



5.1.3 3.3 nm PAH absorption at high temperature 

We compare here 3.3 |i.m observational data with in- 
frared absorption spectra of heated gas-phase PAHs with 
internal temperatures of ~500 K. Data for a signifl- 
cant num ber of small to medium-siz e d PAHs are avail- 
able f rom Semmler. Yang &i Crawford! (|l99ll ). IJoblin et all 



(|1994 I1995I) and from the NIST Standard Refe rence 
Database INIST MS Data Center / S. E. SteinI (|201ll ). In 
Figure 6 infrared absorption spectra of four 4-ring PAH 
compounds are shown: the compact PAH pyrene (CieHio) 
with no bay hydrogens, and three 4-ring isomers (C18H12) 
- benz[a]anthracene, chrysene and triphenylene which have 
2, 4 and 6 bay hydrogens, respectively. The spectroscopic 
data used in th e figure are from the NIST Standard Refer - 
ence Database (|NIST MS Data Center / S. E. SteinI I2OIII ') 
and were recorded at ~500 K. Similar recordings, with three 
molecules i n common and one additional isome r, are given in 
Figure 7 of lSemmler. Yang fc CrawfordI (|l99ll V The spectra 
for a given molecule are an average over many C-H stretch 
vibrational modes, including hot bands, and are remarkably 
simple in appearance. A particularly striking aspect of Fig- 
ure 6 is the shift in band centre to shorter wavelength with 
increase in the number of bay hydrogens on the PAH struc- 
ture. 
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Figure 6. Absorption spectra in the 3.3 ]im region from 
INIST MS Data Center / S. E. SteTnl l|201lh for (top to bottom) 
pyrene (lOnb), benz [a] anthracene (lOnb, 2b), chrysene (8nb, 4b) 
and triphenylene (6nb, 6b). As per convention the hydrogens are 
not shown but those at bay sites are indicated by a star. 



In their study of the absorption sp e ctra o f 33 gas- 
phase PAHs, ISemmle~ Yang fc CrawfordI (|l99lh reported 
the measurement of band wavenumbers, the most significant 
of which for this work are given in Table 2. We have chosen 
to separate their reported wavenumbers into two columns 
labelled Di and 1^2, which have mean values of 3061 and 
3084 cm~^ (3.267 and 3.243 iim), respectively. These values 
represent the summed contributions from many vibrational 
modes of each individual molecule as well as being aver- 
aged over the molecules listed in Table 2. They fall to the 
short wavelength side of the deduced 'astronom ical' com- 
ponents of 3.28 and 3. 30 um dSong et al.ll2007l). However, 
the experimental data of lSemmler. Yang fc CrawfordI (|l99ll ) 
were recorded at ^500 K which is a lower 'temperature' 
than is expected for UV/vis pumping followed by IR cas- 
cade emission in an astronomical context for this size of 
molecule, and so a small shift to longer wavelength is ex- 
pected. In the next Section a model taking into consider- 
ation the emission process of PAH molecules is presented. 
A significant result is that the difference of 0.024 \im in the 
average bay/non-bay laboratory data is close to the wave- 
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length separation of ~0.02 [im obtained for the two compo- 
nents (3.28 and 3.30 [im] fro m analysis of the astronomical 
data (|Song et al.ll2003l , [2007h . It therefore appears that the 
laboratory data fall largely into two frequency categories 
which, we suggest here, follow the degree of bay/non-bay 
character. This has much in common with the idea of an 
infrared 'group frequency' where it is the chemical nature 
of a group that determines the vibrational frequency rather 
than it being specific to an individual molecule. Accord- 
ingly the lower (nb) and higher (b) C-H frequencies should 
hold for most PAH molecules. Although the discussion in 
this Section has focussed on the bay/non-bay spectroscopic 
characteristics, it is important to note that the formation 
of bay-type sites arises naturally during growth of PAHs, 
for example in the growth from the symmetrical compact 
(non-bay) molecule pyrene into benzo[e]pyrene. Hence PAH 
growth and the emergence of bay-type hydrogens are linked. 
In addition the similar spatial distributions of the 3.28 (im 
and 3.40 ^im band carriers (Section l4.1|l and as discussed in 
ISong et al.l (|2007l ) suggests that the 3.28 vim emission may 
also arise from addition of a chemical group (in this case 
C4H2) to the edge of a PAH structure which is consistent 
with formation of a 'bay' configuration. 



6 MODELLING OF THE C-H BAND 

In order to test the bay /non-bay hypothesis, the photo- 
physics of PAH molecules needs to be taken into account. 
After the absorption of a UV photon, the energy of the 
molecule is redistributed in the vibrational modes through 
Intramolecular Vibrational Redistribution (IVR), and the 
molec ule cools by emitting in its IR active m odes (see for ex- 
ample lAUamandola. Tielens fc Barkeij|l989l ). The approach 
used in modelling the 3.3 [im feature consists of two parts: 

a) the vibrational infrared transition frequencies and tran- 
sition strengths of a molecule are calculated using DFT and 

b) these data are used as input to compute the emission 
spectrum. 

6.1 DFT calculations 

PAH vibrational transition frequencies and strengths were 
calculated using the B3LYP functional with a 6-310* basis 
set within the Gaussian software package (|Frisch et al.l 
I2004D . This basis set was selected in preference to the 
more commonly used 4-31G basis as it provides better 
agreement between calculated an d experimental intensitie s 
for the C-H stretching modes (|Pathak fc Rastogil l2007h . 
Comparison of the computed DFT spectra with available 
experimental absorption spectra shows that the DFT 
method overestimates both the absolute C-H frequencies 
(as is well known) and also the frequency difference between 
bay and non-bay C-H stretches. The discrepancies are 
larger when B3LYP/4-31G is used. It is important to 
recognise that there is a significant redshift in the observed 
wavelength for high internal temperatures; for example 
at ~500 K the peak wavelength (Table 2) for pyrene is 
3.28 vim whereas in emission following U V laser-desorption 
excitation the emission peak is at 3.30 \xm (|Schlemmer et ahl 
1 19941 ) in good agreement with astrophysical observation. 
To account for the differences between the DFT and 
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experimental absolute and relative {i.e. bay c.f. non-bay) 
frequencies, two slightly different scaling factors were used 
and applied to all computed C-H stretching modes: 0.948 
for the 2800 - 3210 cm"^ range (non-bay) and 0.945 for 
cj > 3210 cm^^ (bay). The first factor (non-bay Hs) was 
determined by taking the experimental UV las er-desorbed 
experimental result of ISchlemmer et af] (|l994l ) for pyrene 
with the band position-tem p eratur e law for pyrene as 
determined by Ijoblin et al.l (|l995h : a temperature of 
1300 K for 'astrophysical pyrene' is inferred from the 
observation, which corresponds to the absorption of a 
6.35 eV photon in the UV-pumped modelling. The sec- 
ond scaling factor (for bay Hs) was calculated using the 
mean separation between experimental bay and non-bay 
modes from Table [S] It is notable that the experimental 
bay c.f. non-bay frequency difference does not change 
significantly with molecular size for the small-to-medium 
sized PAHs of relevance here. Insp e ction of Table 1 in 
iBauschlicher. Peelers fc AUamandolal (|2009l ) for the C-H 
stretch modes in large PAH molecules (Nc ^ 60) containing 
bay regions shows that the splitting with B3LYP/4-31G 
takes a range of higher values between 0.03 and 0.06 [im 
but these large molec ules are not thought to be the carriers 
of the 3.3 yim band (lAU amandola, Ti elens fc Barked Il989l : 
ISchutte, Tielens fc Allamandola,1993i ). 



6.2 The 3.3 band profile 

The transition frequencies and strengths from the DFT cal- 
culations were used as in put to a UV-pumped PA H emis- 
sion model as employed bv lPathak fc Rastogil (|2008l ). In this 
model the rate of UV photon absorption by a PAH molecule 
is given by: 

Jz.2 

where By is the Planck function of HD 44179 with T^a — 
7750 K (jMen'shchikov et allliooj ). u is the UV photon fre- 
quency and the molecular photo-absorption c ross-section, 
oVjis taken from the French-Italian databasJl llMulas et al.l 
l2006t ). The emission spectrum arising from all bay and non- 
bay IR active C-H stretch transitions is then calculated. A 
Lorentzian with FWHM of 25 cm~^ was used, this being the 
mean value deduced fro m the two-component spectroscopic 
fits o f ISong et al.l(|2007l ). Anharmonic effects, for which there 
is little evidence in the form of 'hot bands' in the astronom- 
ical spectra, and rotational broadening, were neglected. 
Figure [7| shows a comparison between astronomical spectra 
of the Red Rectangle at three offsets and calculated emis- 
sion spectra. Pyrene was chosen as representative of a non- 
bay structure (10 nb; CieHio) and perylene (8 nb and 4 b; 
C20H12) as an example of a slightly larger molecule contain- 
ing bay regions (Figure [SJ. In the top panel of Figure [T] the 
on-star 3.30 \ira b and is fitted well by a pyrene-only con- 
tribution. Indeed, ISong et al.1 (|2003h found that the on-star 
feature can be modelled by a single Lorentzian profile. The 
middle panel shows a coadded spectrum of pyrene and pery- 
lene (1:1.5 ratio) compared with the spectrum of the Red 

^ http:/ /astrochemistry.ca. astro. it/database/ 



© 2012 RAS, MNRAS OOO.mfTOl 



8 A. Candian, T.H. Kerr, I.-O. Song, J. McCombie and P.J. Sarre 



c 



0) 

> 
o 



0.; 



0.0 



0.5 



0.0 



0.5 






RR 

model /A 
/ / V 


on stor 




1/ \ 
f \ 
r \ 
' \ 












H 1 1 1 1 1 1 1 1- 

V \ 


9 no" 












-. — . — 1 ^ 




-, 1 1 , , , 1 , - 

//\ 


D.DD 




/ V 

/ V 
1 V 

7 \ 
/ \ 





3.25 3.30 3.35 
Wavelength (/^m) 



3.40 



Figure 7. Comparison between tlie 3.3 |xm band in the Red Rect- 
angle at 0", 2" and 5.66" offset from HD 44179 (solid line) and 
coadded emission spectra of pyrene and perylene (dashed Une) 
in defined ratios (see text). Astr onomical data ar e taken from 
ISong et all | |2003|) (top panel) and lSong et al.l | |2007|) (central and 
bottom panels). The Pf(5 emission line at 3.296 |xm has been re- 
moved from the on-star spectrum, and is responsible for the slight 
deviation from a smooth profile at this wavelength. The band 
modelling is described in Section 6. 



Rectangle at 2.00" offset along the NW axis. The presence 
of a molecule with bay regions (perylene), permits a very 
good fit of the 3.3 \im band in a part of the nebula where 
the strengt h of the 3.28 \im component starts to grow (c./. 
Figure 3 in ISong et al.|[2007l ). The bottom panel shows the 
astronomical spectrum at 5.66" offset compared with the 
coadded emission of pyrene and perylene in a 1:7 ratio (not- 
ing that perylene has 2:1 non-bay/bay hydrogens). The in- 
creased intensity of the 3.28 \m\ component further into the 
nebula can therefore be accounted for in terms of a more sub- 
stantial contribution from bay-containing structures. This 
does not mean that perylene of itself gives rise to most of 
the 5.66" offset feature but rather that it would indicate an 
approximate ratio between non-bay and bay hydrogens of 
the order of 2:1 in the PAHs present at this offset. As the 
common Type 1 (A3, 3) feature is yet more symmetrical a 
more even ratio of non-bay/bay Hs is anticipated. This re- 
sult is consistent with mass spectrometric analyses of PAHs 
in meteorites which r eveal bay-contai n ing P AHs to be com- 
mon (see for example I Callahan et al.l (|2008l ) and references 
therein). Thus the bay /non-bay combination offers a satis- 



factory new explanation for the two components and pro- 
file variation in the 3.3 \im band. It contrasts strongly with 
the proposal that the 3.3 vim band is a superposition prin- 
cipally of anionic and neutral PAHs, the negatively charged 
PAHs being responsible for the longer w avelength part 
(|Bauschlicher. Peelers fc Allamandolall2009l ): this would re- 
quire PAH anions to be the most abundant PAHs close to 
HD 44179 which seems unlikely. Their ionisation energies 
are less than one eV and so these anions would be sus- 
ceptible to electron photodetachment and hence removal of 
the carrier. A larger set of small PAHs is required to test 
the bay/non-bay interpretation more fully and is beyond 
the scope of this paper. However, there must be some ques- 
tion of uniqueness in the absence of additional information. 
It may also be noted that the transition strengths calcu- 
lated using DFT cannot be fully relied on, perhaps partic- 
ularly for subtle changes such as those between bay and 
non-bay environments; for example in the case of triph- 
enylene (6b and 6nb) the experimental data (Figure 6) 
show the strongest absorption to fall in the the bay spec- 
tral region whereas the relatively high-level DFT calcula- 
tions employed here predict bay and non-bay peaks to be 
of roughly equal intensity. Nevertheless, the fact that the 
observed 3.3 \xm profile and its evolution with offset can 
be accounted for based on experimental data, theoretical 
calculations and emission modelling leads to the conclusion 
that there is no need to invoke positively and/or negatively 
charged PAHs to explain the 3.3 \im f e ature as dis cussed by 
iBauschlicher. Peelers fc AUamandoial (|2008l . |2009| ) . 



7 OTHER RED RECTANGLE EMISSION 

The Red Rectangle is remarkable for its variety of spectral 
features in both absorption and emission across a range of 
wavelengths. Unidentified optical emission bands and Ex- 
tended Red Emission (ERE) appear prominently along the 
whiskers of the nebula. However, the carriers of the rela- 
tively new ly discovered 'blue luminescence' are centred on 
HD 44179 |viih. Witt fc Gordon|[2004l . l2005l ). This emission 
has been attributed to the fiuorescence of optically pumped 
gas-phase PAHs including the small neutral molecules an- 
thracene and pyrene. It might therefore be expected that 
the UIR emission from neutral PAHs and blue luminescence 
would share the same spatial distribution but this is not the 
case. The blue luminescence is strongest near 3934 A and, 
while centrally peaked, it is clearly extended E - W along the 
direc tion of the circumbinary disk (see Figure 1 of lViih et al.l 
l2006h : this suggests a chemical or physical association be- 
tween blue luminescence and the disk. It is notable that the 
spatial distribution at 4050 A is, in contrast, highly symmet- 
rical about the central star and is likely due predominantly 
to scattered starlight. 

At ~ 5"S offset there is reasonable agreement be- 
tween the 3.3 \im and blue luminescence spatial distributions 
(|Viih. Witt fc Gordon 2005 ). but in the inner 3" region 
they are markedly different. It therefore appears that the 
carriers of blue luminescence and 3.3 [im emission in the in- 
ner region are not the same entities. This peculiarity could 
possibly be resolved if the blue luminescence were to arise 
from PAHs such as anthracene and pyrene adsorbed on cool 
grains in the disk (~ E-W), whereas the 3.3 jxm emission 
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were due to gas-phase PAHs entrained in the ~ N-S outflow. 
A recent near-IR polarimetric study of the inner part of the 
Red Rectangle has revealed scattering at J which traces the 
circumbinary disk material and which appear s also to follow 
the s patial distribution of blue luminescence (|Gledhill et al.l 
|2009| '1. A measurement of the polarisation properties of blue 
luminescence would be of particular interest in assessing 
whether the blue luminescence carriers are free gas-phase 
molecules or molecules adsorbed on grains. 



of iBauschlicher. Peeters fc AUamandolal (|2009h for the Red 
Rectangle where the 3.3 ^i.m feature is attributed to a mix- 
ture of anionic, neutral and (marginal) cationic PAHs. 

In this paper we have described new IFU observations of 
the spatial distribution of PAHs that give rise to the 3.3 ^.m 
emission band and a proposed two-component interpreta- 
tion. The images provide support for the proposal, based 
originally on spectroscopic long-slit data, that there are two 
types of C-H carrier that we suggest can be accounted for 
in terms of non-bay and bay hydrogen sites. 



8 DISCUSSION AND CONCLUSIONS 



In earlier work (|Song et al.ll2003l . l2007^ it was proposed that 
the variation in profile of the 3.3 vim PAH band with offset 
in the Red Rectangle could be interpreted in terms of the 
size distribution of neutral PAHs. Following the additional 
IFU observations reported here, and based on our analysis 
of additional laboratory data, consideration of the spectro- 
scopic properties of bay/non-bay structures and emission 
modelling, we propose it is the emergence of bay-type struc- 
tures of growing PAHs that underlies the observed profile 
change. The 3.3 [un emission band of the Red Rectangle may 
then be interpreted in terms of compact non-bay PAHs dom- 
inating on-star with emission centred at 3.30 \im, with larger 
PAH molecules which include bay-type structures (and a 
typical wavelength of 3.28 ii.m) increasing in abundance with 
offset from HD 44179. The emergence of bay-type structures 
and the issue of PAH size are connected as the growth of 
PAHs must pass through stages where stoichiometric C4H2 
entities are added to a PAH. Long slit spectra in the 10- 
15 ^im region in different positions of the nebula are needed 
to fully understand the variation of molecular edges, i. e. bay 
and non-bay regions, and to compare with the different spa- 
tial distribution of the two components of the 3.3 [im band. 

We have considered other possible explanations for 
two 3.3 \im emission band components including the ef- 
fects of hydrogenation and methylation which are ex- 
pected to increase with offset from the star HD 44179, 
and give rise to the 3.4 [im emission band. This band 
has a similar spatial distributio n to that o f the 3 .28 i-i-m 
component - see Figure 2 and ISong et af] (|2003l . |2003). 
However, we have not found a systematic effect of ei- 
ther hydrogenation or methylation of PAH molecules on 
the 3.3 [im aromatic C-H frequency in available l abora- 
tory and theoretical data JWagner. KinTfc Saykall^ 2000l : 
iPauzat. Talbi fc EmngeJll999l : iPauzat fc Ellingedl200lf r 



Bay-type sites in a more general astronomical context 
can arise by addition of further rings as a PAH structure 
grows or alternatively through damage to the edge of a com- 
pact PAH. We favour the former interpretation for the in- 
ner part of the Red Rectangle given the good correlation in 
growth of the 3.28 (im and 3.4 pim features with offset, as 
the 3.4 \im band is generally accepted as arising from addi- 
tion of hydrogen and/or methyl groups. In other astronom- 
ical environments such as H 11 regions lik e Orion H2S1 the 
observed very strong 3.28 ii.m component (|Song et al.ll2007l ') 
could arise because of damage to the edge structure. In the 
interpretation put forward here, the widely observed A3. 3 
(Type 1) 3.3 \un profile refiects a superposition of transi- 
tions of small neutral PAHs with bay (3.28 \im) and non-bay 
(3.30 (i.m) edges; it differs considerably from the conclusion 
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